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We explore stability of gravitating Q-balls with potential V^{(p) = ^-0 2 — A0 4 + -^- via catastrophe 

theory, as an extension of our previous work on Q-balls with potential V-j,{(f>) = ^-<p 2 — ^(/) 3 + \(fi 4 . 
In flat spacetime Q-balls with V4 in the thick- wall limit are unstable and there is a minimum charge 
Qmin, where Q-balls with Q < Q m in are nonexistent. If we take self-gravity into account, on the 
other hand, there exist stable Q-balls with arbitrarily small charge, no matter how weak gravity is. 
That is, gravity saves Q-balls with small charge. We also show how stability of Q-balls changes as 
gravity becomes strong. 
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I. INTRODUCTION 

Q-balls [lj, a kind of non-topological solitons Q, ap- 
pear in a large family of field theories with global U(l) 
(or more) symmetry, and could play important roles in 
cosmology. For example, the Minimal Supersymmetric 
Standard Model may contain baryonic Q-balls, which 
could be responsible for baryon asymmetry Q and dark 
matter Q. 

Because Q-balls are typically supposed to be mi- 
croscopic objects, their self-gravity is usually ignored. 
Therefore, stability of Q-balls has been intensively stud- 
ied in flat spacetime 0t3| • Q-balls in arbitrary dimension 
Q and spinning Q-balls 0, [k} have also been studied. 

If Q-balls are so large or so massive, on the other hand, 
their size becomes astronomical and their gravitational 
effects are remarkable [Hi [HI- For example, it has been 
shown [lj|] that the size of Q-balls is bounded above due 
to gravity. There are analogous objects which are anal- 
ogous to gravitating Q-balls: boson stars [li|. While 
Q-balls exist even in flat spacetime, boson stars are sup- 
ported by gravity and nonexistent in flat spacetime. Al- 
though a difference in theory between Q-balls and boson 
stars is solely the potential parameters, investigations of 
their properties have been carried out separately so far. 

In our previous paper [l4[ , to obtain a unified picture of 
Q-balls and boson stars, we made an analysis of gravitat- 
ing Q-balls and boson stars via catastrophe theory (l5j . 
In Ref. 1J| we chose a potential for Q-balls 



V 3 {((>) := — 2 -^(/> 3 +A(/> 4 , with m 2 , fj,, A > 0, (1.1) 
because in the limit of fi — > this approaches a typical 



potential for boson stars, 



Vbs(4>) ■= + M\ with 



A>0. (1.2) 
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As a result, we found that Q-balls and boson stars expose 
a similar phase relation between a charge and a total 
Hamiltonian energy. (See, cusp structures in Figs. 1(a) 
and 12(a) in [3.) 

In this paper we extend our analysis via catastrophe 
theory to a potential 

W) : = — tf-^+jfi with ™ 2 ' A ' M > °> (°) 

which we call V4 Model [l6j]. We choose this potential 
because previous work on Q-balls in flat spacetime @, 0] 
showed stability of Q-balls with V3 Model (|1.1[> and V4 
Model (|1.3I) are quite different. We are interested in how 
gravitating Q-balls properties depend on potentials and 
what universal properties are. 

This paper is organized as follows. In Sec. II, we de- 
rive equilibrium field equations. In Sec. Ill, we show nu- 
merical results of equilibrium Q-balls and discuss their 
stability. In Sec. IV, we discuss why thick-wall solutions 
become stable against the naive expectation that gravity 
is not effective for Q-balls with small charge. In Sec. V, 
we devote to concluding remarks. 



II. ANALYSIS METHOD OF EQUILIBRIUM 
Q-BALLS 

A. Equilibrium field equations 

We begin with the action 
S = / d^x^J—g 
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(2.1) 
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where <p 



is an S0(2)-symmetric scalar field 



and <f> := \J <p ■ cb — y/4>\+ §\- We assume a spherically 
symmetric and static spacetimc, 



ds 2 



-a 2 (r)dt 2 +A 2 (r)dr 2 +r 2 (d9 2 +sm 2 0dip 2 ). (2.2) 



For the scalar field, we assume that it has a spherically 
symmetric and stationary form, 



{4>ii4>2) = 4>(r) (cos wt, sin ut). 
Then the field equations become 
rA 3 



(2.3) 
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(2.6) 



where ' := d/dr. To obtain Q-ball solutions in curved 
spacetime, we should solve (|2.4|) - ()2.6p with boundary 
conditions, 

A(0) = A(oo) = a(oo) = f , 

A'(0) = a'(0) = 0'(O) = 0(oo) = 0. (2.7) 

We also restrict our solutions to monotonically decreasing 
0(r). Due to the symmetry, there is a conserved charge 
called Q-ball charge, 



Ar 2 (f> 2 

where / := 4tt / dr. 



as 



(2.8) 

We suppose V4 Model (|1 ,3[) . Rescaling the quantities 

t := XMt, f := XMr, 
Va rh 2 72 



m 
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\/AM' 
, k := GAM 2 , (2.9) 



the field equations (|2.4[) - (|2.6[) with the potential (|1.3p are 
rewritten as 

* + >-i)-w(& + g + v 4 ), 

(2.10) 



(2.11) 



» + £_4^ + ^y#_A"«. (2.12) 



B. Stability analysis method via catastrophe theory 

In our previous paper [l4[ , we discussed how we apply 
catastrophe theory to the Q-ball and boson star systems. 
Here, we summarize our method. An essential point is 
to choose behavior variable(s) , control parameter's) and 
a potential in the Q-ball system appropriately. 

We use the Hamiltonian energy £ as a potential be- 
cause 8E/84>\q = 8E/8g^ v — 0, reproduces the equilib- 
rium field equations (I2.4|) - (12.6[) . The Hamiltonian energy 
E was calculated as [1 4j 



r-»oo 2GA 



Ms 
2 ' 



(2.13) 



where Ms is the Schwarzschild mass. We also normalize 
E and Q as 



E 



M 



Q := QX. 



(2.14) 



Because the charge Q and the model parameters m 2 
and k specify the system environment, they should be 
regarded as control parameters. To discuss a behavior 
variable we consider a one-parameter family of perturbed 
field configurations <f>x(/) near the equilibrium solution 
4>(r). Because dE[<j> x ]/dx — (8E/8cf> x )d(f> x /dx = when 
4> x is an equilibrium solution, a; is a behavior variable. 

According to Thorn's theorem, if the system has two 
control parameters, there is essentially one behavior vari- 
able; if the system has three control parameters, there are 
one or two behavior variables. Because the present Q- 
ball system contains (Q,?ti 2 ,k), we speculate that each 
has two behavior variables, ui 2 and 0(0). However, be- 
cause stability structure of equilibrium solutions in three- 
parameter space (Q,m 2 ,K) is very complicated and our 
interest is how gravitational effects change the stabil- 
ity structure, in what follows, we discuss the stability 
structure of equilibrium solutions in two-parameter space 
(Q, k) under fixed rh 2 . 

Our method of analyzing the stability of Q-balls is as 
follows. 

• Fix the value of rh 2 . 

• Solve the field equations (|2.4[) - (I2.6[) with the bound- 
ary condition (|2.7p numerically to obtain equilib- 
rium solutions 4>(r) for various values of uj and k. 

• Calculate Q for each solution to obtain the equilib- 
rium space M = {(x, Q, k)}. We denote the equa- 
tion that determines M by f(x,Q,K) = 0. 
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Find folding points where dQ/dx = or dn/dx = 
0, in M, which are identical to the stability-change 
points, S = {(x, Q, k) \ df/dx = 0, / = 0}. 



Calculate the energy E by (I2.13|) for equilibrium so- 
lutions around a certain point in E to find whether 
the point is a local maximum or a local minimum. 
Then we find the stability structure for the whole 
M. 



III. EQUILIBRIUM SOLUTIONS AND THEIR 
STABILITY 



In preparation for discussing gravitating Q-balls, we 
review their equilibrium solutions and stability in flat 
spacetime (k = 0). The scalar field equation (|2.12j) re- 
duces to 



77/ 2 7/ -2 7 . dv i 
b =—-<j> =- . 



(3-1) 



This is equivalent to the field equation for a single static 
scalar field with the potential := V4— w 2 0/2. Equilib- 
rium solutions satisfying boundary conditions (|2.7I) exist 
if min(K,) < 14(0) and (PV^ / 'd<f> 2 (0) > 0, which is equiv- 
alent to 
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0<e <r 



(3.2) 



where e := \Jfin 2 



The two limits e 2 



\ and e40 

correspond to the thin- wall limit and the thick- wall limit, 
respectively. 

It is usually assumed that the potential has an abso- 
lute minimum at <f> = 0. If V(0) is a local minimum but 
the absolute minimum is located at <f> 7^ 0, true vacuum 
bubbles with charge (Q-bubbles) may appear. The con- 
dition for Q-bubbles is rh 2 < 0.5. Therefore, stability 
structure falls into two classes, rh 2 < 0.5 and rh 2 > 0.5 
0: 

• m 2 > 0.5: For each m 2 , there is a nonzero minimum 
charge, Qmim below which equilibrium solutions do 
not exist. For Q > Qmin, stable and unstable solu- 
tions coexist. 

• m 2 < 0.5: For each rh 2 , there is a maximum charge, 
Omax, as well as a minimum charge, <2 m i n , where 
one stable and two unstable solutions coexist for 

Qmin ^ Q ^ Qm&x- For Q <C Qm'm Or Q 7> Q ma x, 

there is one unstable solution. 

To discuss gravitational effects later, it is useful to es- 
timate the central value 0(0) in flat spacetime. Because 
V u = at spacial infinity, its order of magnitude is esti- 
mated as a solution of V u = (0(0) 7^ 0). For V4 with 
the thick- wall condition e < 1, we obtain 



Thus, 0(0) - e. 

It was shown [lj| that in V3 Model properties of grav- 
itating Q-balls also depend on whether rh 2 > 0.5 or 
rh 2 < 0.5. In the following analysis, therefore, we choose 
rh 2 = 0.6 and 0.3 typically. Other cases are not qual- 
itatively different from these cases. For our numerical 
calculation, we use the Bulirsch-Stoer method based on 
the double precision FORTRAN program. 
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FIG. 1: Behavior of the metric A(r) and <j>(r) for the solutions 
cj 2 ~ 0.34 in (a) and (b), respectively. We find that the scalar 
field is concentrated near the origin as shown in (b). This 
tendency becomes stronger as gravity is stronger. Thus, A(r) 
varies near the origin compared with that in the thick-wall 
solutions. 



A. Gravitating Q-balls for m 2 > 0.5 

In this subsection we fix m 2 = 0.6. First, we present 
typical solutions in Figs. [1] and [2] we choose k = 0, 
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(a) 



FIG. 2: Behavior of the metric A(r) := y/g r r and cj>(r) for the 
thick- wall solutions Co 2 ~ 0.595 in (a) and (b), respectively. 



0.03, and 0.1 and show the metric A(r) in (a) and the 
scalar field amplitude <j>(r) in (b). In Figs. Q] we put 
Cj 2 ~ 0.34. We find that as k becomes larger, or gravity 
is stronger, |^4(r) 2 — 1| becomes up to order one, and the 
Q-ball size becomes smaller by self-gravity. As we shall 
discuss below, the solutions with k = 0, 0.03 in Figs. [Hare 
stable, while the solution with k = 0.1 is unstable. That 
is, strong gravity destabilizes or kills some of the solutions 
which would be existent and stable without gravity. 

Figs. [2] show the solutions with w 2 ~ 0.595. Because 
e 2 = 0.05 <C 1, these are thick-wall solutions. We find 
an interesting feature in (b): the difference between <f> 
with k — 0.03 and n = 0.1 are small, but they are quite 
different from </>(r) with k = 0. This indicates that the 
configuration of <j>(r) for gravitating Q-balls does not ap- 
proach that for k = if we take the limit of k — > 0. In 
the next section we shall discuss the reason for this. 

In this way we calculate equilibrium solutions numer- 
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FIG. 3: (a) Q-E and (b) Q-u 2 relations for rh 2 = 0.6 with k = 
0, 0.03 and 0.1. For the flat case K = 0, it has been found that 
solutions with solid (dotted) lines are stable (unstable) [?J. 
We extend these interpretations for the gravitating case. 



ically for various w 2 and show Q-E and Q-u; 2 relations 
in Figs. [3J We can obtain stability of the solutions using 
catastrophe theory as follows. 

• When there are multiple values of E for a given set 
of the control parameters (to 2 , k, Q), by energetics 
the solution with the lower value of E should be 
stable. 

• Once the stability for a given set of the parameters 
(m 2 ,K,Q) is found, the stability for all sets of pa- 
rameters which are reached continuously from that 
set without crossing turning points (i.e., stability- 
change points E) is the same. 

• Stability changes across E. 

• Spiral structure in the Q-uj 2 plane should be consid- 
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FIG. 5: Stability interpretation via catastrophe theory for 
m 2 = 0.6 for n = 0.03. 



FIG. 4: Stability interpretation via catastrophe theory for 
rh 2 = 0.6 for the flat case. 



ered exceptionally. We interpret that all solutions 
are unstable there. 

As a result, we can conclude that solid and dashed lines 
correspond to stable and unstable solutions, respectively. 

To illustrate this energetic or catastrophic argument 
more clearly, we give a sketch of the potential function E 



near the equilibrium solutions in Figs. [4] and [5] Figs. [4] 
shows the flat case (k = 0): the solid (dotted) line in (a) 
and the points B, a± and b\ correspond to those in (b), 
where candidates of behavior variables (f>(0) and Co 2 are 
shown, (c) is a schematic picture of the potential function 
E in terms of the behavior variable uj 2 . The point B, the 
cusp in (a), corresponds to the inflection point in (c); we 
understand why there is no solution with lower Q. The 
point a\ on the solid line and the point b\ on the dotted 
line in (ii) correspond to the potential minimum and the 
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FIG. 6: Structures of the equilibrium spaces, M = 
{(w 2 , k, Q)}, and their catastrophe map, x(M), into the con- 
trol planes, C = {(ft, Q)}, for in 2 = 0.6. Blue lines and red 
lines in M represent stable and unstable solutions, respec- 
tively. In the regions denoted by Si, SiU (i = 1,2) and N 
on C, there are one stable solution, i stable solution(s) and 
one or more unstable solution(s), and no equilibrium solution, 
respectively, for fixed (ft, Q). 



maximum in (c), respectively. 

As another example for catastrophic interpretation, 
we discuss stability of the solutions for k = 0.03, using 
Figs. [5] A complicated structure appears in the enlarge- 
ment in (a). In the Q-range between A and B there are 
triple values of E for fixed Q. In this case the poten- 
tial function should be given by (ii) in (c). This means 
that two stable solutions coexist for fixed Q in this range. 
As a result, we can conclude that there are stable grav- 
itating Q-balls which approach Q — > in the thick-wall 
limit (ui 2 — > 0.6). Fig. [5] (b) tells us that in the unstable 
sequence right the point C there is no one-to-one cor- 
respondence between ui and the solutions while there 
is one-to-one correspondence between </>(0) and the solu- 
tions. In this range, therefore, <fr(0) is more appropriate 



as a behavior variable than Co 2 as shown in (iii) and (iv) 
in (c). 

(a) 
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FIG. 7: Q-E relations for m 2 = 0.3 with k = and 0.05 in 
(a) and with k = 0.1 and 0.12 in (b), respectively. As in the 
case for m 2 = 0.6, gravitating cases have sequences from A 
to the origin, which are considered to be stable, written by 
solid lines. For ft = 0.12, the sequence from B-C-D separates 
into two parts. Each branch has sequences of cusp structures 
(around Q ~ 8 and 48 in (b)). 

Fig. [6] shows the structures of the equilibrium spaces, 
M = {(Cj 2 , k,Q)}, and their catastrophe map, x(M), into 
the control planes, C — {(k,Q)}, for rh 2 = 0.6. In the 
regions denoted by Si, S^U (i = 1, 2) and N on C, there 
are one stable solution, i stable solution(s) and one or 
more unstable solution(s), and no equilibrium solution, 
respectively, for fixed (k,Q). The points A, B and C 
in Figs. [3] are marked by circles, triangles and squares, 
respectively. For example, for n — 0.1, which is the case 
shown in Figs. [31 we can confirm that only a stable solu- 
tion exists below Q ~ 4 (the point D) which is denoted 
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FIG. 8: Q-u) 2 relation for the same solutions as in Figs. [7] The 
points from A to D and the dotted (solid) lines correspond 
to those in Figs. [7] From this correspondence, we find that 
flat Q-balls approach Q — > oo in the thick-wall limit ui 2 — > 0.3 
while gravitating Q-balls approach Q — > in this limit as in 
the case for in 2 — 0.6. 

by Si in Fig. [6l One stable solution and one or more 
unstable solution(s) exist in the region from Q ~ 4 to 
Q ~ 8 (the point A) which is denoted by SiU. 

Main characteristics of the equilibrium solutions in 
Figs. [3] and [6] are summarized as follows. 

• If re = 0, there is a minimum charge, Q m ; n , denoted 
by B on the re = line in Fig. [3] The equilibrium 
solutions in the thick-wall limit e 2 — > are unsta- 
ble, as indicated by the dotted lines. 

• If n ^ 0, no matter how small K is, the equilibrium 
solutions in the thick-wall limit e 2 —> are stable 
and Q — > 0. These stable solutions correspond to 
the solid lines from Q = to A in Fig. [3] We can 
interpret that gravity saves thick-wall Q-balls. 

• If k ^ 0, the maximum charge, Q max , emerges in 
the thin- wall range. This extreme solution is de- 
noted by C on the k = 0.03 line in Fig. [3] That is, 
gravity kills Q-balls with large charge. 

• If re ^ 0, spiral trajectories appear in the Q-Cj 2 
plane. 

• If re ~ 0.1, the two extremal solutions B and C 
merge and disappear, and accordingly the stability 
sequence between them disappear, too. 

The second result is remarkable. In Sec. IV, we investi- 
gate why these discontinuous changes occur at re = 0. 
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FIG. 9: Stability interpretation via catastrophe theory for 
m 2 = 0.3 and re = 0.05. 



B. Gravitating Q-balls for m 2 < 0.5 

In this subsection we fix m 2 = 0.3. We show Q-E 
relations in Figs. [7] and Q-Cj 2 relations in Fig. [3] In the 
same method as in Sec. Ill A, we can determine stability 
of the equilibrium solutions: solid lines and dashed lines 
correspond to stable and unstable solutions, respectively. 

In the case of re = 0, for example, there are two cusp 
structures as shown in Figs.[7J(a). Only solutions in the 
narrow range between B and C are stable. As another 
example, for re = 0.05 we illustrate catastrophic inter- 
pretation in Figs. IHl In the Q-range between A and B 
there are quadruple values of E for fixed Q. In this case 
the potential function is given by (ii) in (c). As in the 
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FIG. 10: Structures of the equilibrium spaces, M = 
{((1> 2 , k, Q)}, and their catastrophe map, \(M), into the con- 
trol planes, C = {(k, Q)}, for m 2 = 0.3. Blue lines and red 
lines in M represent stable and unstable solutions, respec- 
tively. In the regions denoted by SiU (i = 1, 2) and N on C, 
there are i stable solution(s) and one or more unstable so- 
lution^), and no equilibrium solution, respectively, for fixed 

(«,<?)• 



solution(s). 

Main characteristics of the equilibrium solutions in 
Figs. [JJ [5] and [TUJ are summarized as follows. 

• If k = 0, there is a maximum charge for stable so- 
lutions, Qmax, denoted by C, as well as a minimum 
charge, Q m im denoted by B, on the k = line in 
Figs. [JJ [5] The equilibrium solutions in the thick- 
wall limit e 2 — > are unstable, as indicated by the 
dotted lines. 

• If ft ^ 0, no matter how small k is, the equilibrium 
solutions in the thick-wall limit e 2 — > are stable 
and Q — > 0. These stable solutions correspond to 
the solid lines from Q = to A in Figs. [7J We can 
interpret that gravity saves thick-wall Q-balls. 

• As k increases, the maximum charge, <5 max , in- 
creases until k = K cr ; t ~ 0.1. 

• The solution sequence for fixed k splits into two 
when k = K cl -it. In each sequence spiral trajectories 
appear in the Q-ui 2 plane. 



IV. THICK- WALL LIMIT 

It is not surprising that properties of gravitating Q- 
balls change gradually as k increases. It seems strange, 
however, properties of gravitating Q-balls in the limit of 
k — > differs completely from that of flat Q-balls (n = 0), 
as show in Figs. 151 and ITUl Here we discuss the reason for 
this. 

We consider the case of weak gravity (k <C 1) and 
thick- wall (e 2 1). Since the gravity is weak, we can 
express the metric functions as 

a 2 = l + h(r), A 2 = l + /(r), (h < 1, / « 1). (4.1) 

Up to first order in h and /, we can rewrite the scalar 
field equation (j2.6l) as 



case of to 2 = 0.6 and k = 0.03 in Fig. [5l two stable so- 
lutions coexist for fixed Q in this range, and there are 
stable gravitating Q-balls which approach Q — > in the 
thick-wall limit {Hi 2 — > 0.3). 

Figure [TU] shows the structures of the equilibrium 
spaces, M = {(uj 2 ,k,Q)}, and their catastrophe map, 
x(M), into the control planes, C = {(k, Q)}, for m 2 
I)..'! x(M) shows that in the regions denoted by S^U 
(i = 1, 2) and N on C, there are i stable solution(s) and 
one or more unstable solution(s), and no equilibrium so- 
lution, respectively, for fixed (k,Q). The points A, B 
and C in Figs. [JJ are marked by circles, triangles and 
squares, respectively. For example, for n = 0.05, if we fix 
Q between ~ 12 (the point B) and ~ 18 (the point A), 
there are two stable solutions and one (or more) unstable 



+ ~ - £) $' = (l+/)[( e 2 +^ 2 )0 2 -40 3 +60 5 ] . 

If we fix e 2 > and take the limit of n —> (i.e., 
h, f —> 0), Eq. (|4.2[) reduces to the field equation in 
flat spacetime, (|3.1j) . However, if we take the limit of 
e 2 — > as well, the situation becomes complicated. For 
any small n, if we take so small e 2 that e 2 <§; hQ 2 , which 
means e 2 kw 4 as we shall show below, the first or- 
der term in h, hu 2 (j) 2 , dominates the zeroth order term, 
e 2 cf) 2 . That is, in the thick-wall limit of e 2 — > 0, the scalar 
field equation with infinitesimally small k can be different 
from that with k = 0. 

The above argument is based on the hypothesis that 
h does not approach zero as fast as e 2 when we take the 
limit of e —> 0. Otherwise, the inequality e 2 < hui 2 would 
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be wrong. To complete this argument, we shall estimate 
the order of magnitude of h by assuming 

4>{f) ~ e < 1 for r <-, (4.3) 

which is valid in flat spacetime. From the Einstein equa- 
tions, we find 

(4.4) 

where i runs spacial components. If we take the weak 
field approximation (|4.ip and the thick-wall approxima- 
tion e 2 <C 1, we obtain 

(r 2 ti)' ~ 87r K f 2 m 2 2 . (4.5) 

With the boundary condition h'(0) = and the approx- 
imation (14.31) . we can integrate (|4.5I) as 

8 I 
h! ~ — 7TKTO 2 e 2 f, for f < — . (4-6) 
3 e 

With the boundary condition h(oo) = and the approx- 
imation (14.31) , we can integrate (|4.6I) as 

/i ~ ^7TKTO 2 (e 2 f 2 - I ), for f<-. (4.7) 

This means h ~ Km 2 , which is independent of e. There- 
fore, we can conclude e 2 <C /iw 2 in (|4.2[) in the thick-wall 
limit with fixed k. 

If the assumption (14.3[) is not valid, the configuration of 
</>(r) is quite different from that for flat spacetime. This 
case also means that gravitating Q-balls are completely 
different from those in flat case even if the gravity is very 
weak. 

From the above argument, we can understand why 
gravity saves thick-wall Q-balls. This is not surprising 
because similar phenomenon occurs in the case of boson 
stars: boson stars with Vbs do not exist in flat case while 
they exist even if the gravity is very weak. 

V. CONCLUSION AND DISCUSSION 

We have analyzed stability of gravitating Q-balls with 
V\ via catastrophe theory. Our results are summarized 



as follows. 

Although our original concern was massive Q-balls 
with astronomical size, we have found an unexpected re- 
sult that the weak gravity changes properties of thick- 
wall Q-balls. In flat spacetime Q-balls in the thick- wall 
limit are unstable and there is a minimum charge <5 m i n , 
where Q-balls with Q < Q m [ n are nonexistent. If we take 
self-gravity into account, on the other hand, there exist 
stable Q-balls with arbitrarily small charge, no matter 
how weak gravity is. That is, gravity saves Q-balls with 
small charge. 

This result indicates that gravitational effects may be 
important for other models, such as Q-balls in supersym- 
metric extensions of the Standard Model. For example, 
gravity may allow for a new branch of solutions in some 
parameter range where equilibrium solutions are nonex- 
istent in the absence of gravity. 

We have also shown how stability of Q-balls changes 
as gravity becomes strong. For example, if m 2 > 0, the 
maximum charge, Q ma x, decreases as gravity becomes 
strong, while there is no maximum charge in flat space- 
time. That is, gravity kills thin-wall Q-balls with large 
charge. 

In the case of strong gravity, only Q-balls with small 
charge exist, and instability solutions make spiral trajec- 
tories in the Q-Cj 2 plane. These properties are common 
to Q-balls with V3 potential and boson stars with Vbs- 
While Q-balls and boson stars have been studied sepa- 
rately so far, our result suggest that there is universal 
nature of gravity, which may be important to discuss Q- 
balls with astronomical size or boson stars. 
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